Bound on the time variation of the fine structure constant 
driven by quintessence 
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The bound on the time variation of the fine structure constant (a) driven by the dynamics of 
quintessence scalar field which is coupled to electromagnetism is discussed using phenomenological 
quintessential models constrained by SNIa and CMB observations. We find that those models 
allowing early quintessence give the largest variation Ace at the decoupling epoch. Furthermore, the 
fifth force experiments imply that Act/a is less than about 0.1%. 
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PACS numbers: 98.80.Cq, 98.80.Es, 06.20.Jr 

The idea of time-varying fundamental physical con- 
stants such as gauge coupling constants was proposed 
long time ago [l| . Recently it has been claimed that the 
results of a search for time variability of the fine structure 
constant (a), using absorption systems in the spectra of 
distant quasars, yield 5-a evidence for a smaller a in the 
past: Aa/a = -0.543 ± 0.116 x 10~ 5 over the redshift 
range 0.2 < z < 3.7 Q. In addition, there exist terrestrial 
and cosmological constraints. An analysis of the isotropic 
abundances in the Oklo natural uranium fission reac- 
tor, active about 1.8 x 10 9 years ago (corresponding to 
z ~ 0.14), suggests -0.9 x 10~ 7 < Aa/a < 1.2 x 10" 7 Q. 
A bound Aa/a < 3 x 10~ 7 has been obtained from the 
analysis of the Re/Os ratio in meteorites (z ~ 0.45) 
While the data of cosmic microwave background (CMB) 
anisotropics are consistent with a being smaller by a 
few percent at the decoupling epoch (z ~ 1100) in a 
flat cold dark matter model with a cosmological con- 
stant (ACDM) @, the recent observations made by the 
Wilkinson Microwave Anisotropy Probe (WMAP) have 
provided a bound -0.06 < Aa/a < 0.02 at 95% CL @. 
At much higher redshifts, big bang nucleosynthesis con- 
siderations place bounds |Aa/a| < 10~ 2 at z ~ 10 9 0. 

The recent astrophysical and cosmological observations 
such as type la supernovae (SNe) and CMB anisotropics 
concordantly prevail the ACDM model @- The current 
data, however, are consistent with a somewhat broader 
diversity of such "dark energy" as long as its equation 
of state (EOS) approaches that of the cosmological con- 
stant at a recent epoch. A dynamically evolving scalar 
field </> called "quintessence" (Q) is probably the most 
popular scenario so far to accommodate the dark energy 
component. Many Q models have been proposed with 
various physically motivated effective potentials V(<f>) for 
the scalar field Q. Several attempts have been made to 
test different Q models ^(J. Nevertheless, it proves to be 
premature at this stage to perform a meaningful data fit- 
ting to a particular Q model, or to differentiate between 
the variations. The reconstruction of V(4>) would likely 
require next-generation observations. As such, model in- 
dependent approaches which simply involve parametriz- 
ing the dark energy EOS have been proposed to comply 
with the SNe and CMB observational data 



It is of particular interest to study the observational ef- 
fects of direct interaction of <f> to ordinary matter if there 
is any. For example, imposing an approximate global 
symmetry would allow a coupling of <f> to electromag- 
netism, which would lead to rotation of polarized light 
from distant radio sources, temporal evolution of a |14j . 
and generation of primordial magnetic fields [l5]. Re- 
cently, there have been many studies on the time- varying 
a in the context of quintessential cosmology b y in vok- 
ing non-renormalizable </>-photon couplings [m, Il7l [l8 j . 
However, most of the studies are based on certain Q 
models |l7| . In this paper, we will discuss the observa- 
tional constraints to the evolution of <f> through its EOS 
P<p = W4>P4> in the parametrized Q models [H], 0, ^| 
that in turn give rise to a model independent bound on 
the Q-driven time variation of a at the decoupling epoch. 

The </>-photon coupling that wc consider here is 
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where F^„ = d^A v — d v A^, M p is the reduced Planck 
mass (87rG)~ 1 / 2 , and k is a coupling constant. The fifth- 
force experiments limit 0] 



k < 1Q- 
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Therefore, by defining 8 = <f}/M p , the change of a from 
the present time to a time t is simply given by 



^=0 - 9(t). 
Ka 
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Consider a flat universe in which the total density param- 
eter of the Universe today is represented by fio = + 
n° r +n° = 1 with 9? m h 2 = 0.135 and n° r h 2 = 4.152x 10~ 5 , 
where the present Hubble constant is parametrized as 
Hq = 100/ikms _1 Mpc _1 . Assuming a spatially homoge- 
neous 4> field, the evolution of the cosmic background is 
governed by 
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where a is the scale factor and the conformal time 
is defined as rj — H J dta^ 1 ^). Here we have used 
(d(j)/dt) 2 = (1 + w^p^ and V(<p) = (1 - w <f> )p <f> /2, and 
rescaled the energy density of the ith component as 
Qi = pi j (MpH ) 2 ■ Accordingly, the dimensionless Hub- 
ble parameter is given by 

n 2 = (h/h q ) 2 = n° m a - 3 + n° a - 4 + n 2 . (6) 

Therefore, given the prescribed EOS w$(r\) which is a 
function of conformal time, the problem is reduced to 
solving a set of first-order coupled ordinary differential 
equations l@J and (JSJ) with the initial conditions set at 
the present time and the time variable running backward. 
Let us introduce a useful quantity that is the f^-weighted 
average [2^ 

rvo / rno \ - 1 

(wj,) = / n^w^drj x ( / fyWdr,) , (7) 
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where 770 and ?7dcc are the conformal time today and at 
decoupling respectively. It is known that 7^ af- 
ter decoupling would result in a time- varying Newtonian 
potential which produces large-scale CMB anisotropy 
through the integrated Sachs- Wolfe (ISW) effect. It was 
shown 20, 21] that as far as the CMB anisotropy power 
spectrum is concerned, the time- averaged (w^) can be 
well approximated by an effective constant value for w$ 
as long as the Q component is negligible at decoupling. 

Physically, — 1 < w<f, < 1, where the former equal- 
ity holds for a pure vacuum state. Lately some progress 
has been made in constraining the behavior of <p from 
observational data. A combined large scale structure, 
SNe, and CMB analysis has set an upper limit on Q 
models with a constant < —0.7 |2J, and the 
recent WMAP CMB data gives a stronger limit <w$ < 
—0.78 23]. Furthermore, the SNe data and measure- 
ments of the position of the acoustic peaks in the CMB 
anisotropy spectrum have been used to put a constraint 
on the present w° < —0.96 |24|. The apparent brightness 
of the farthest SN observed to date, SN 1997ff at redshift 
z ~ 1.7, is consistent with that expected in the deceler- 
ating phase of the flat ACDM model with tt A ~ 0.7 |25|. 
implying w^, = — 1 for z < 1.7. 

Let us first work on the simplest case with a constant 
W0 = —0.8 to obtain an order-of-magnitude estimate on 
Aa/na. By solving 



d9/drj = ay (I + w^qj,, (8) 

we find that Aa/ra decreases quickly with redshift to a 
constant value of —0.82 at z ~ 20. Next we study the 
varying a in two phenomenological Q models which are 
based on the maximum likelihood analysis of the current 
SNe and CMB observational data. The first model has a 
monotonic EOS w^a) = wq + (w m — wq)(1 — a) (where 
W0 — * w m as a — > 0) that has been shown to be able to 
describe the late-time (f> evolution for a wide class of Q 
models |ll|. Using wq = —0.99 and w m — —0.31, which 



are within 2a from the best-fit parameters to the SNe and 
CMB data [2l|. we find that Aa/na decreases to —1.63 
at decoupling. The second model has an EOS with four 
parameters W = (wo,w m , a c , A) |l2^ : 

(1 + e a '=/ A ) (1 - eM- a ^ A ) 
w<t>(a) = w + (w m -w ) ^ - e(ac _ a)/A ) ^ _ gl/A ) , (9) 

where a c signifies the transition epoch and A represents 
the transition rate from wq to w m . We have worked out 
two cases: a slow EOS with W = (-1, -0.3,0.1,0.8) and 
a rapid EOS with W = (-1,0,0.4,0.1), both of which 
are taken from the la values of the best-fit models to the 
WMAP and SNe data [27j]. We find that Aa/na decrease 
to —1.56 and —3.82 at decoupling for the slow and rapid 
EOS respectively. In fact, the slow EOS and the first 
model are very similar, in which n^/D, m decreases to 0.01 
at z ~ 30 and becomes negligible in the dark ages. On 
the contrary, the rapid EOS which has a larger change 
in a allows early quintessence with J7^,/J7 m ~ 0.1 for 
1100 > z > 10. This can be understood by considering a 
constant for which the problem has a simple solution. 
Equation Q gives cx a _3 ( 1+ "'*) and Sl ( j > /Q, m oc a -3 ™*. 
Thus, from Eq. JSJ we have dO/d-q oc a~( 1+3 ™'*)/ 2 and 
therefore the change of a will grow if w,/, > —1/3. Also, 
the condition for early quintessence is simply > 0. 
From the results of the above slow and rapid EOS's where 
w m = —0.3 and w m — respectively, we see that the 
change of a in models with > > —1/3 is smaller 
than that in models with early quintessence. Therefore, 
early quintessence is a crucial condition for a large tem- 
poral change of a. We thus see that within viable models 
the varying a strongly depends on the Q component in 
the dark ages. 

To further study the Q dependence and to find out 
the maximum change of | Aa/ na\, we turn to the generic 
quintessence (GQ) model [L|, in which a piecewise con- 
stant EOS subject to observational constraints is intro- 
duced. The GQ model starts with determining the form 
for Wcf, at low redshifts. We choose w<f, ~ — 1 for z < 2 to 
satisfy the above-mentioned SNe constraints. Then 
should be increased to 1 after z = 2 to have a maximum 
change of a based on Eq. (JSJ). However, to avoid an un- 
acceptably large ISW effect, must drop to —1 at a 
certain redshift. This results in the square-wave EOS in 
the GQ model. A GQ model is shown in Fig.^ where the 
width of the square wave is determined by (w^) = —0.8 
to saturate the WMAP upper limit on a constant [23[ . 
In this case, Aa/na reaches a minimum value of —0.4. In 
fact, a square-wave or pulse-shaped EOS is anticipated 
and quite general in the class of Q models using pseudo 
Nambu-Goldstone boson fields incorporated with quan- 
tum corrections to the cosine potential 28] . However, one 
can still shift the whole square wave to higher redshifts 
while fixing (w^) = —0.8 as to generate the required ISW 
effect. We find that the shift increases the Q component 
at the decoupling epoch and affects the location of the 
CMB acoustic peaks. This can be understood as follows. 
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FIG. 1: Aa/na in the GQ model with a square-wave EOS 
located at low redshift. 

The tightly coupled baryon-photon plasma experi- 
enced a serial acoustic oscillation just before the recom- 
bination epoch. The acoustic scale (the angular momen- 
tum scale of the acoustic oscillation) sets the location of 
the peaks in the CMB anisotropy power spectrum [2^| . 
and is characterized by 

I A = nd*/h s = tt(t] - Tydcc)/ / c s dr) , (10) 

Jo 

where represents the comoving distance to the de- 
coupling epoch and h s denotes the sound horizon at 
decoupling, both of which are affected in the pres- 
ence of quintessence. The sound speed c s in the pre- 
recombination plasma is given by |29| 

c s = l/x/3(l + R) with 

4p 7 ^2.725if J 1 + 2 ' v ; 

where the baryon-photon momentum density ratio R sets 
the baryon loading to the acoustic oscillation of CMB. 
Hence, h s at decoupling can be determined by the differ- 
ential equation 

dh s /dri = c s , (12) 

coupled with the background evolution equations and 
©. Using n° b h 2 = 0.0224, h = 0.71, T° = 2.725JT, and 
Zdcc = 1089, we can calculate various acoustic scales by 
shifting the square wave of toward the decoupling 
epoch while keeping (w$) fixed. Figure [3 plots the re- 
sults against the Q field energy density at decoupling 
^ cc with contours (u^) = -0.7, -0.8, -0.9, and -1. 
The last one is simply the ACDM model. All GQ models 
lying in the region above the curve of {w^) = —0.8 and 
within the WMAP bounds are consistent with the cur- 
rent CMB data. Thus, the constant w$ = —0.8 model 



FIG. 2: Acoustic scales of the GQ models are plotted as a 
function of the quintessence density at decoupling while keep- 
ing (Wtf,) fixed. The model equivalent to the cosmological 
constant case has Ia — 301 and is denoted by an arrow. The 
two horizontal lines signify the 1-a upper and lower bounds 
permitted by the WMAP data 
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FIG. 3: Ao/kq in the extreme GQ model with ({vuq) , O.^) 
= (—0.8, 0.26). The decoupling epoch is marked as the ver- 
tical line at Zd cc = 1089. 



and that in Fig. ^ are in fact disfavored as far as the 
acoustic peak location is concerned. In Fig. [31 we have 
worked out an extreme model with parameters taken at 
the upper right corner of the allowed region in Fig.|5J The 
model shows that the Q component dominates over the 
matter density at high redshifts. This early quintessence 
drives Acv/kcv to a minimum value of —2.77. Note that a 
non-zero £1^ after decoupling would lead to a suppression 
of the growth of matter perturbations on scales smaller 
than the smoothing scale of the quintessence. Following 
Ref. 01 > we have found that the growth function for mat- 
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ter perturbations in the extreme GQ model is still at an 
acceptable level with the COBE data. Furthermore, per- 
turbations in the dark energy component will mildly af- 
fect our results because dark energy with (w^,) = —0.8 is 
already close to a vacuum state |30j. At last, the extreme 
GQ model with f^ oc = 0.26 is consistent with the upper 
bound f^0 < 0.39 during the radiation dominated epoch 
obtained by performing a maximum likelihood analysis 
on the CMB data HJ- 

In conclusion we have investigated the quintessence- 
induced evolution of a allowed by the observational con- 
straints from CMB and SNe using the parametrized 
quintessence models. Although the true equation of 
state, if there is any, may be a complicated function 
of time, they should capture the generic features of 
quintessence evolution for a model-independent study 
of the time- varying a. We have found an extreme GQ 
model and a rapid equation of state that allow early 
quintessence driving a maximum change of \Aa/na\ ~ 
3 — 4 at the decoupling epoch. Hence, the fifth-force 
limit in Eq. J2J implies that |Aa/a| is less than about 
0.1% at decoupling as long as the quintessence is coupled 



to photon. Future CMB data will be able to constrain a 
up to 0.1% level or tighter than this when combined with 
other cosmological measurements |(| . 

So far, we have not considered using the phenomeno- 
logical quintessence models to fit the observational data 
of varying a at low redshifts (z < 3.7). Although it 
is ad hoc and fine-tuned, we find that it is not so diffi- 
cult to modify the equation of state at low redshifts to 
make the change of a consistent with the measurements 
of quasar absorption spectra and the Oklo and Re/Os 
limits. Indeed, using a model-independent approach it 
was shown that these low-redshift measurements can con- 
strain the dark energy equation of state today to sat- 
isfy — 1 < w<p < —0.96 and disfavor late-time changes 
in |18| . This result is consistent with what we have 
assumed for the low-redshift value of in our GQ mod- 
els and the constraint Aa/a = —0.06 ± 0.06 x 10~ 5 re- 
cently obtained based on a new sample of high-redshift 
(0.4 < z < 2.3) quasar absorption line systems |32fl . 
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